Proton- 3 He elastic scattering at low energies and the "A y Puzzle" 
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Abstract. The Kohn variational principle and the hyperspherical harmonic technique are applied to study 
p - 3 He elastic scattering at low energies. Preliminary results obtained using several interaction models are 
reported. The calculations are compared to a recent phase shift analysis performed at the Triangle University 
Nuclear Laboratory and to the available experimental data. Using a three-nucleon interaction derived from chi- 
ral perturbation theory at N 2 LO, we have found a noticeable reduction of the discrepancy observed for the A y 
observable. 



1 Introduction 

The four nucleon system has been object of intense studies 
in recent years. In first place, this system is particularly in- 
teresting as a "theoretical laboratory" to test the accuracy 
of our present knowledge of the nucleon-nucleon (NN) 
and three nucleon (3N) interactions. In particular, the effect 
of (i) the NN P-wave and of (ii) the 3N force are believed to 
be larger than in the A — 2 or 3 systems. Moreover, it is the 
simplest system where the 3N interaction in channels of to- 
tal isospin T = 3/2 can be studied. In second place, there is 
a number of reactions involving four nucleons which are of 
extreme importance for astrophysics, energy production, 
and studies of fundamental symmetries. As an example, 
reactions like d + d — > 4 He + y or p + 3 He — > 4 He + v e + e + 
(the hep process) play important roles in solar models and 
in the theory of big-bang nucleosynthesis. 

Nowadays, the four-nucleon bound state problem can 
be numerically solved with good accuracy. For example, 
in Ref. Q me binding energies and other properties of 
the ff-particle were studied using the AV8' ||2) NN interac- 
tion; several different techniques produced results in very 
close agreement with each other (at the level of less than 
1%). More recently, the same agreement has also been ob- 
tained considering different realistic NN+3N force mod- 
els II3I4I5I6I . 

In recent years, there has also been a rapid advance in 
solving the four nucleon scattering problem with realistic 
Hamiltonians. Accurate calculations of four-body scatter- 
ing observables have been achieved in the framework of 
the Faddeev-Yakubovsky (FY) equations [7 8 9], solved in 
momentum space, where the long-range Coulomb interac- 
tion is treated using the screening-renormalization method 
111011 II . Also solutions of the FY equations in configu- 
ration space 11211311411 and several calculations using the 
resonating group model I15I16I17I18I were reported. Cal- 
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culations of scattering observables using the Green Func- 
tion Monte Carlo method are also underway [ 19 1. 

In this contribution, the four-body scattering problem 
is solved using the Kohn variational method and expanding 
the internal part of the wave function in terms of the hyper- 
spherical harmonic (HH) functions (for a review, see Ref. 
[20 1). Previous applications of this method B21I22I13I were 
limited to consider only local potentials, as the Argonne 
uig (AV18) (231 NN potential. Recently, the HH method 
has been extended to treat also non-local potentials, given 
either in coordinate- or momentum-space [24 25 1. A first 
application of the HH method to compute n — 3 H elas- 
tic observables with non-local potentials was reported in 
Ref. [26 1 . Here, we report the first results obtained for p - 
3 He. 

The potentials used in this paper are the I-N3LO model 
by Entem and Machleidt (27), with cutoff A = 500 MeV, 
and the AVI 8. The I-N3LO potential has been derived us- 
ing an effective field theory approach and the chiral pertur- 
bation theory up to next-to-next-to-next-to-leading order. 
We have also performed calculations by adding to the I- 
N3LO potential a 3N interaction, derived at next-to-next- 
to leading order in Refs. [1281291 (I-N3LO/N-N2LO inter- 
action model). The two free parameters in this N-N2LO 
3N potential have been chosen from the combination that 
reproduces the A — 3,4 binding energies ll29l . The devel- 
opment of a 3N interaction including N3LO contribution is 
still under progress 1 30 1 . Finally, we have also considered 
the Urbana IX (UIX) 3N potential El, which has been 
used together with the AVI 8 potential (AV18/UIX interac- 
tion model). 

The four-body studies performed so far have eviden- 
tiated several discrepancies between the theoretical pre- 
dictions and experimental data. Let us consider first n - 
3 H elastic scattering. Calculations based on NN interac- 
tion models only disagree 1211 121321 161 13171 rather sizably 
with the measured total cross section 1551 . both at zero en- 
ergy and in the "peak" region (E„ ss 3.5 MeV). Such an 
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observable is found to be very sensitive to the NN interac- 
tion model 0. At low energy, the discrepancy is removed 
by including a 3N force fixed to reproduce the triton bind- 
ing energy B21I12I16I . but it remains in the peak region. 
Interestingly, this disagreement is noticeably reduced us- 
ing the I-N3LO/N-N2LO interaction model |26|. A similar 
situation occurs in the differential cross section, but defini- 
tive conclusions in this case are difficult to extract since the 
experimental errors are rather large. 

In this respect, the p - 3 He elastic scattering is more 
interesting since there exist several accurate measurements 
of both the unpolarized cross section 041351361 and of the 
proton analyzing power A y [37 22 36]. The calculations 
performed so far (with a variety of NN interactions, and 
with the AV18/UIX model) have shown a glaring discrep- 
ancy between theory and experiment for A y [ 32 2 21 16136171 . 
This discrepancy is very similar to the well known "A y 
Puzzle" in N-d scattering. This is a fairly old problem, al- 
ready reported about 20 years ago [38 39] in the case of n— 
d and later confirmed also in the p -d case Il40l . The inclu- 
sion of standard models of the 3N force has little effect on 
these observables. To solve this puzzle, speculations about 
the deficiency of the NN potentials in 3 Pj waves (where 
the spectroscopic notation 2s+l Lj has been adopted) have 
been advanced. The situation of other p - 3 He observables 
(the 3 He analyzing power Aq v and some spin correlation 
observables as A yy , A xx , etc.) is less clear due to the lack 
(until recently) of equally accurate measurements. 

However, recently 0411421 at the Triangle University 
National Lab. (TUNL) there has been a new set of accu- 
rate measurements (at E p = 1.60, 2.25, 4 and 5.54 MeV) 
of various spin correlation coefficients, which has allowed 
a phase-shift analysis (PSA). The aim of this paper is to 
compare the results of the theoretical calculations to these 
new data. Moreover, we want to study the effect of includ- 
ing the N-N2LO 3N force in A y . 

This paper is organized as follows. In Section|2] a brief 
description of the method is reported. In Section[3] a com- 
parison between HH and FY calculations is shown. We 
have performed this comparison for n - 3 H scattering using 
the I-N3LO potential for incident neutron energy E„ = 3.5 
MeV. In Section |4] the theoretical calculations are com- 
pared with the available experimental data. The conclu- 
sions will be given in Section 



region where the particles are close to each other and their 
mutual interactions are strong. On the other hand, W^f de- 
scribes the relative motion of the two clusters in the asymp- 
totic regions, where the p - 3 He interaction is negligible 
(except for the long-range Coulomb interaction). In the 
asymptotic region the wave functions iPj^L reduces to , 
which must therefore be the appropriate asymptotic solu- 
tion of the Schrodinger equation. can be decomposed 
as a linear combination of the following functions 



4 

1=1 

x \fi(yi) ±i 
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where yi is the distance between the proton (particle t) and 
3 He (particles ijk), q is the magnitude of the relative mo- 
mentum between the two clusters, s/ the spin state of par- 
ticle I, and 03 is the 3 He wave function. The total kinetic 
energy T cm in the center of mass (cm.) system and the 
proton kinetic energy E p in the laboratory system are 
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^ p ^ L cm. 



(3) 



where p = (3/4)Mjy is the reduced mass of the system 
(Mn is the nucleon mass). Moreover, Fi and Gl are the 
regular and irregular Coulomb function, respectively, with 
rj = 2pe 2 /q, and D^+i is a normalizing factor defined to be 
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The function fi_(y) — \X~ exp( _ y%)] 2i+1 in Eq. has been 
introduced to regularize Gl at small y, and fiiy) — * 1 as y 
is large, thus not affecting the asymptotic behavior of fff 3 . 
Note that for large values of qyi, 
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where <r L is the Coulomb phase shift. Therefore, Q^ s (Q LS ) 
describe the asymptotic outgoing (ingoing) p- 3 He relative 
motion. Finally, 



2 The HH Technique for Scattering States 

In the following, we consider a /?- 3 He scattering state with 
total angular momentum quantum number JJ Z , and parity 
n (the dependence on the wave function and other quanti- 
ties on JJ,n will be understood in the following). The wave 
function f^ 3 describing the two particles with incoming 
relative orbital angular momentum L and channel spin S 
(S = 0, 1) can be written as 



7\ 
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where the parameters SLs^s'(q) are the S -matrix elements 
which determine phase-shifts and (for coupled channels) 
mixing parameters at the energy T cm _. Of course, the sum 
over L' and S ' is over all values compatible with the given 
/ and parity n. In particular, the sum over L' is limited to 
include either even or odd values such that (-1) L = n. 

The "core"wave function W^f has been here expanded 
using the HH basis. For four equal mass particles, a suit- 
able choice of the Jacobi vectors is 



where the part Yf? vanishes in the limit of large inter- 
cluster separations, and hence describes the system in the 
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where p specifies a given permutation corresponding to the 
order i, j, k and / of the particles. By definition, the permu- 
tation p = 1 is chosen to correspond to the order 1, 2, 3 
and 4. 

For a given choice of the Jacobi vectors, the hyper- 
spherical coordinates are given by the so-called hyperra- 
dius p, defined by 



\ P + x lp + x i P 



(independent of p) , (8) 



and by a set of angular variables which in the Zernike and 
Brinkman |43 44| representation are (i) the polar angles 
Xj p = (0j p , (pip) of each Jacobi vector, and (ii) the two addi- 
tional "hyperspherical" angles ip2 P and ip 3p defined as 
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where Xj p is the modulus of the Jacobi vector Xj. The set 
of angular variables x\ p , X2 P , x 3p , <p2 P , and fop is denoted 
hereafter as Q p . The expression of a generic HH function 
is 

«,K,AM (Q ) = 



X},, 
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where P" n ' b are Jacobi polynomials and the coefficients A/fj^ 3 
normalization factors. The quantity K = l\ + I2 + (3 + 
2(«2 + «3) is the so-called grand angular quantum number. 
The HH functions are the eigenfunctions of the hyperangu- 
lar part of the kinetic energy operator. Another important 
property of the HH functions is that p K 'H^'^ Lj ni n (& P ) 
are homogeneous polynomials of the particle coordinates 
of degree K. 

A set of antisymmetrical hyperangular-spin-isospin states 
of grand angular quantum number K, total orbital angu- 
lar momentum A, total spin E, and total isospin T (for the 
given values of total angular momentum J and parity n) 
can be constructed as follows: 

^ = Z*f CT (U,*.0, (ID 

where the sum is over the 12 even permutations p = i, j, k, I, 
and 
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Here, < Hf l j , 2 % M ^ (Q P ) is the HH state defined in Eq. QH, 
and Si (tj) denotes the spin (isospin) function of particle i. 
The total orbital angular momentum A of the HH function 
is coupled to the total spin E to give the total angular mo- 
mentum JJ Z , whereas n = (-l) f[+f2+e3 . The quantum num- 
ber T specifies the total isospin of the state. The integer 
index p labels the possible choices of hyperangular, spin 
and isospin quantum numbers, namely 

p = {A, h, €3, h, » 2 , n 3 ,S a ,S b , T a , T b } , (13) 

compatibles with the given values of K, A, E, T, J and 71. 
Another important classification of the states is to group 
them in "channels": states belonging to the same chan- 
nel have the same values of angular (\,l2, (3, La, A, spin 
S a ,S b,E, isospin T a , 7^, T quantum numbers but different 
values of «2> "3- 

Each state ^ AZT entering the expansion of the four 
nucleon wave function must be antisymmetric under the 
exchange of any pair of particles. Consequently, it is nec- 
essary to consider states such that 



0^ T (i,j;k;D^-0^ T (j,i;k;l). 
which is fulfilled when the condition 
£ 3 +S a + T a =odd, 



(14) 



(15) 



is satisfied. 

The number Mkaxt of antisymmetrical functions >f^ ALT 
having given values of K, A, E, and T but different com- 
binations of quantum numbers p (see Eq.([T3l)) is in gen- 
eral very large. In addition to the degeneracy of the HH 
basis, the four spins (isospins) can be coupled in differ- 
ent ways to E (T). However, many of the states Y^^ 7 , 
p = 1, . . . , Mkaxt are linearly dependent between them- 
selves. In the expansion of Yff it is necessary to include 
only the subset of linearly independent states, whose num- 
ber is fortunately noticeably smaller than the correspond- 
ing value of Mkaxt- 

The internal part of the wave function can be finally 
written as 



KAXT 



(16) 



KAXT p. 



where the sum is restricted only to the linearly independent 
states. We have found convenient to expand the "hyperra- 
dial" functions uj^j, (p) in a complete set of functions, 
namely 

M-l 

u KAZTp(P) = C KAXTp,m 9miP) , (17) 
m=0 



and we have chosen 
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(m + 8)! 
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where Lf\yp) are Laguerre polynomials [45] and y is a 
parameter to be variationally optimized. 
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The main problem is the computation of the matrix el- 
ements of the Hamiltonian. This task is considerably sim- 
plified in two steps. First, by using the following transfor- 
mation 

^ T (i,j;k;D = Yja™F(p)&K, Ai:T (l,2;3;4) , (19) 

where the coefficients a^^f T (p) have been obtained using 
the techniques described in Ref. ll46l . Second, by "project- 
ing" the asymptotic states over a complete set of angular- 
spin-isospin states, constructed in terms of the Jacobi vec- 
tors Xj corresponding to the particle order 1 , 2, 3, 4, as fol- 
lows: 



(y f3 (* 3 )Oi*2)s 2 ). (^ 2 (£>)* 3 ). 



2 F^ ± (x u x 2 , x 3 )y a (x u x 2 , x 3 ) , (20) 

a 

where 

(21) 

r- 

and a = {€i,€ 2 , t 3 , 73, Ji,S2, T2, T 3 , T). Note that due 
to the antisymmetry of the wave function, we must have 
{3 + S2 + T2 — odd. This "partial wave expansion" is per- 
formed to include all states a such that l t < £ m - dx . The 
functions F^ s± can be obtained very accurately by direct 
integration 



F^ ± (x l ,x 2 ,x 3 ) 



I 



dx\dx2dx 3 



y a {x\,X2,h) 



Q ± 



(22) 

using Gauss quadrature techniques. 

Finally, the wave function using the transformation ( TT9b 
and the partial wave expansion given above can be rewrit- 
ten as 



1+3 



2 (xi ,x 2 , x 3 )J/ ff (ii , x 2 , X3) 



(23) 



where T is a combination of Jacobi polynomials of the hy- 
perangles, functions g m {p), and functions F a coming from 
the asymptotic parts. Then, the matrix elements of a two- 
body potential can be evaluated as following. First of all 

<^" 3 I Tj V ^f>\ ¥ M ,S > = 6(¥^ 3 \V(l,2W^') , (24) 

due to the asymmetry of the wave function. The matrix 
element of V(l, 2) (note that x 3 = r2-rj) is easily obtained 
using the decomposition given in Eq. ( 1231 . For example, 
for a non-local potential 

<^ + s 3 |V(i,2)|<^'> = 

= J" d 3 X\d 3 x 2 d 3 x 3 d 3 x' 3 ^F^ 3 (xi, X2, x 3 )j 



Table 1. Phase-shift and mixing angle parameters for n - 3 H 
elastic scattering at incident neutron energy E„ = 3.5 MeV cal- 
culated using the I-N3LO potential. The values reported in the 
columns labeled HH have been obtained using the HH expansion 
and the Kohn variational principle, whereas those reported in the 
columns labeled FY by solving the FY equations |8 |. 
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The calculation of the above integral is performed in two 
steps. First, the spin-isospin-angular matrix elements 



dx\dx.2dx-idx' 3 }/ a (x\, x 2 , x 3 ) 

xV(cc3,:e 3 ) y a '(x\,X2,x' 3 ) 



are computed analytically, and, second, the integration over 
the moduli of the Jacobi vectors, 

J" dx\dx2dx 3 dx' 3 x\x^x\x' 3 {j :r {f (xi,X2,x 3 )J 

xy i£:3f (x3 ' x 3) r^ s \x u x 2 ,x^) ,(27) 

can easily obtained by using Gauss quadrature methods, 
being a 4-dimensional integral. The accuracy of the ma- 
trix elements so calculated depends mainly on £ m - ix , the 
maximum value of the orbital angular momentum used to 
truncate the expansion ( 1201 . Values £ m . dx = 5 or 6 have 
been found appropriate to obtain a sufficient numerical ac- 
curacy. 

The S -matrix elements Sls,l'S' (p) an d coefficients m 
occurring in the expansion of are determined by mak- 
ing the functional 

[SLS,L<s-(q)] = SLs,L<s>(q) - j. ( F&f ' \H - E\ (28) 



and 

KAZT(i,m 



stationary with respect to variations in the Sls,us 
(Kohn variational principle). By applying this principle, a 
linear set of equations for Slsjjs' an d cI £ai:th m * s obtained. 
The linear system is solved using the Lanczos algorithm. 

The main difficulty of the application of the HH tech- 
nique is the slow convergence of the basis with respect to 
the grand angular quantum number K. This problem has 
been overcome by dividing the HH basis in classes. More 
details of this method can be found in Ref. ll36l. 



xV(x 3 ,x' 3 )'FF,1'(x 1 ,X2,x' 3 ) 



(25) 



3 Comparison between HH and FY results 

The calculated phase-shift and mixing angle parameters 
for n - 3 H elastic scattering at E„ = 3.5 MeV using the 



19 International IUPAP Conference on Few-Body Problems in Physics 



o 
-10 

-20 

Tit 
o 

a -30 

| -40 
ft 

-50 
-60 
-70 



1 I 1 I 1 I 1 

TUNL 2009, single-energy 
1-N3LO 

I-N3LO/N-N2LO 

AV18 

AV18/U1X 




1 



2 3 4 5 6 
E [MeV] 



-10 - 
-20 - 
-30 - 
-40 - 
-50 - 
-60 - 
-7oL 




2 3 4 
E [MeV] 



Fig. hp- 3 He S-wave phase shifts for four potential models. 
The results of the PSA performed at TUNL have been also re- 
ported t42l . 

I-N3LO potential are reported in Table Q] The values re- 
ported in the columns labeled HH have been obtained us- 
ing the HH expansion and the Kohn variational principle, 
whereas those reported in the columns labeled FY by solv- 
ing the FY equations JS). As can be seen, there is a good 
overall agreement between the results of the two calcula- 
tions. Similar comparisons for p - 3 He are currently under- 
way. 



4 Results 

In this Section, we present preliminary calculations of elas- 
tic p - 3 He observables. Final results with complete tests 
of convergence with respect to K and { ma - dx will be reported 
in a forthcoming paper. 

In the energy range considered here (E p < 6 MeV), 
the various p - 3 He observables are dominated by S-wave 
and P-wave phase-shifts (D-wave phase shifts give only a 
marginal contribution, and more peripheral phase shifts are 
negligible). Let us first concentrate on the comparison of 
calculated S-wave and P-wave phase shifts with those ob- 
tained by the recent PSA @2J, shown in Figs.QJg] The [ S 
and 3 S i phase shifts are reported in Fig.Q] The results ob- 
tained including NN interactions only slightly overpredict 
(in absolute value) the PSA values. Including the 3N force, 
the calculated phase shifts agree very well with the PSA 
values (for both AV18/UIX and I-N3LO/N-N2LO mod- 
els). In fact, the p - 3 He interaction in S-wave is repul- 
sive, being dominated by the Pauli repulsion, and the cor- 
responding phase shifts are generally well reproduced by 
an interaction model giving the correct value of the 3 He 
binding energy (and radius). 

In Fig. |2] the comparison is shown for the J" — 1~ 
phase shifts. The PSA values for the l P\ phase shift have 
large error bars and all potential models give results con- 
sistent with them. In general, we observe that the inclusion 
of the 3N interaction reduces such a phase shift. For the 
3 Pi phase shift, we observe that the theoretical results ob- 
tained with AVI 8 and AV18/UIX models disagree with the 
PSA values, in particular for E p > 4 MeV. Moreover, the 
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Fig. 2. As in Fig.[TJbut for the i P l and 3 P t phase shifts. 



inclusion of the UIX 3N force diminishes the phase shift, 
enlarging the discrepancy. On the contrary, the results ob- 
tained with I-N3LO and I-N3LO/N-N2LO are very close 
to the data. The inclusion of such a 3N force slightly in- 
creases the 3 Pi phase shift. 

A similar behavior can be observed for the 3 f2 and 
3 fo phase shifts, reported in Fig. [3] The 3 P2 phase shift is 
well reproduced by I-N3LO and I-N3LO/N-N2LO models, 
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Fig. 3. As in Fig.[T]but for the 3 Pi and 3 Pq phase shifts. 
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Fig. 4. As in Fig.[TJbut for the 1 and 1 + mixing parameters. 
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Fig. 5. p - 3 He differential cross sections calculated with the I-N3LO (blue dashed line), the I-N3LO/N-N2LO (blue solid line), and 
the AV18/UIX (thin green solid line) interaction models for three different incident proton energies. The experimental data are from 
Refs. I34I35I36I . 




Fig. 6. p - 3 He A y observable calculated with the I-N3LO (blue dashed line), the I-N3LO/N-N2LO (blue solid line), and the AV18/UIX 
(thin green solid line) interaction models for three different incident proton energies. The experimental data are from Refs. [37 22 36]. 
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while with AVI 8 and AV18/UIX we observe an underpre- 
diction. For the 3 Po phase shift we observe that the models 
including NN interaction only overpredict the PSA values. 
With the inclusion of the 3N forces, the results come close 
to the PSA values. This behavior is not unexpected, as it 
was already known that this phase shift "scales" with the 
3 He binding energy [7]. 

Finally, in Fig. [4] the comparison is performed for the 
1" and 1 + mixing parameters. For s(l~), the PSA results 
have quite large errors, however it is significant that only 
the I-N3LO/N-N2LO values are in the error bar. For s(l + ), 
the results obtained with all the four models are more less 
consistent with the (largely scattered) PSA values. 

In conclusion, with AVI 8, the S-wave and most of the 
P-wave phase shifts are at variance with the results of the 
PSA. Including the UIX 3N force, only the S-wave shifts 
are adjusted. With I-N3LO, the S-wave phase shifts are 
overpredicted (in absolute value), but the P-wave phase 
shifts are much better reproduced. Including the N-N2LO 
3N interaction, the S-wave shifts become well reproduced, 
while the agreement with the P-wave shifts is not spoiled. 
In this case also the mixing parameters are well repro- 
duced. We observe therefore an overall good agreement 
between the results obtained with the I-N3LO/N-N2LO 
model and the PSA values. 

Let us now compare the results of the calculations di- 
rectly with the available experimental data. The potential 
models considered hereafter are the I-N3LO, I-N3LO/N- 
N2LO, and AV18/UIX. The calculated p - 3 He differen- 
tial cross sections a energies E p = 2.25, 4, and 5.54 MeV 
are reported in Fig. [5] and compared with the experimental 
data of Refs. |34 35 36 1. As can be seen, there is a good 
agreement between the theoretical calculations and exper- 
imental data. In fact, this observable is sensitive to small 
changes of the phase shifts only in the "interference" re- 
gion around 9 — 30 deg. We note only a slight deviation of 
the AV18/UTX results from the data in such a region. 

More interesting is the situation for the proton vector 
analyzing power A y . Here, we observe a larger sensitiv- 
ity to the employed interaction model. The calculations 
performed using AV18/UIX largely underpredict the ex- 
perimental points, a fact already observed before B22I361 . 
Using I-N3LO, no significative changes are observed. A 
sizable improvement is found by adopting the I-N3LO/N- 
N2LO model, as it was expected from the discussion re- 
garding the comparison with the PSA phase shifts. To be 



Table 2. Ratio A* eo 7A^ xp at the peak for the four potential mod- 
els and three values of proton energy E p . 



theor /^expt 



at 



more quantitative, let us consider the ratio A h ,/A 
the peak, reported in Tabled For the AV18, AV18/UIX, 
and I-N3LO models, there is approximately a 30% un- 
derprediction of the peak height. The underprediction is 
reduced approximately by a factor 2 for the I-N3LO/N- 
N2LO model. 

Finally, in Figs. Q and [8] we show two further polariza- 
tion observables, the 3 He analyzing power A^o and the spin 
correlation coefficient A m . These observables (and other 
spin correlation coefficients) are not very sensitive to the 
interaction model. We observe however, that the I-N3LO/N- 
N2LO model produces slightly better agreement with data. 
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Fig. 7. p - 3 He observable A QlJ calculated with the I-N3LO (blue 
dashed line), the I-N3LO/N-N2LO (blue solid line), and the 
AV18/UIX (thin green solid line) interaction models for three dif- 
ferent incident proton energies. The experimental data are from 
Refs. (371411 . 
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Fig. 8. p - 3 He observable A m calculated with the I-N3LO (blue 
dashed line), the I-N3LO/N-N2LO (blue solid line), and the 
AV18/UIX (thin green solid line) interaction models for three dif- 
ferent incident proton energies. The experimental data are from 
Refs. (3714T1 . 



5 Conclusions 

In this paper preliminary results for p - 3 He elastic scatter- 
ing obtained using the HH function expansion have been 
reported. We have considered four interaction models. The 
first two are based on phenomenological NN interactions 
(AV18 and AV18/UIX). The other two are derived from 
a chiral effective theory (I-N3LO and I-N3LO/N-N2LO). 
Both NN interactions (AVI 8 and I-N3LO) reproduce the 
NN scattering data with a x 1 P er datum very close to 1 . 
The UIX 3N interaction has been fixed in order to repro- 
duce the trinucleon binding energies, whereas the N-N2LO 
by fitting simultaneously both A — 3, 4 binding energies. 
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Fig. 9. p - d observable A tJ at E p = 3 MeV calculated with the I- 
N3LO (blue dashed line), I-N3LO/N-N2LO (blue solid line), and 
AV 1 8/UIX (thin green solid line) interaction models for E p = 3 
MeV. The experimental data are from Ref. (47l . 



We have compared the results obtained with the four 
interaction models with the available experimental data and 
a very recent PSA performed at TUNL [42] . We have found 
that the phase shifts obtained with I-N3LO/N-N2LO are 
consistent with those derived from the PSA. The direct 
comparison of the I-N3LO/N-N2LO results with the ex- 
perimental data has shown that there are still some dis- 
crepancies, but the Ay problem is noticeably reduced. In 
fact, we observe that now the discrepancy is of the order of 
15% at the peak, much less than what observed with other 
interaction models. 

In this regard, it is interesting to note that in the N - d 
case, the use of the I-N3LO/N-N2LO model does not give 
such a reduction of the N - d "A y Puzzle", as can be seen 
for example in Fig. [9] We observe in fact only a small 
increase of the A y computed with the I-N3LO/N-N2LO 
model with respect to that calculated, for example, using 
the AVI 8/UIX model. Note that also for p - d, the calcu- 
lated A y with I-N3LO and AVI 8/UIX are very similar. In 
this case the ratio A* eor /A^ p at the peak is approximately 
0.81 (0.74) for I-N3LO/N-N2LO (AVI 8/UIX). Investiga- 
tions to understand this behavior are in progress. 

Acknowledgments. We would like to thank TV Daniels 
and T.B. Clegg for providing us with their experimental 
data prior to publication, and A. Deltuva and A. Fonseca 
for useful discussions. 



References 

1. H. Kamada et al, Phys. Rev. C 64, (2001) 044001 

2. B. S. Pudliner et al, Phys. Rev. C 56, (1997) 1720 

3. A. Nogga et al, Phys. Rev. C 67, (2003) 034004 

4. R. Lazauskas and J. Carbonell, Phys. Rev. C 70, 
(2004) 044002 



5. M. Viviani, A. Kievsky, and S. Rosati, Phys. Rev. C 
71, (2005) 024006 

6. R.B. Wiringa et al, Phys. Rev. C 62, (2000) 014001 

7. A. Deltuva and A. C. Fonseca, Phys. Rev. C 75, (2007) 
014005 

8. A. Deltuva and A. C. Fonseca, Phys. Rev. Lett. 98, 

(2007) 162502 

9. A. Deltuva and A. C. Fonseca, Phys. Rev. C 76, (2007) 
021001 

10. E. O. Alt, W. Sandhas, and H. Ziegelmann, Phys. Rev. 
C 17, (1978) 1981; ibid.. 21, (1980) 1733 

11. A. Deltuva, A. C. Fonseca, and PU. Sauer, Phys. Rev. 
C 71, (2005) 054005; ibid., 72, (2005) 054004 

12. F. Cieselski and J. Carbonell, Phys. Rev. C 58, (1998) 
58; F. Cieselski, J. Carbonell, and C. Gignoux, Phys. 
lett. B447, (1999) 199 

13. R. Lazauskas et al, Phys. Rev. C 71, (2005) 034004 

14. R. Lazauskas, Phys. Rev. C 79, (2009) 054007 

15. H. M. Hofmann and G. M. Hale, Nucl. Phys. A613, 
(1997)69 

16. B. Pfitzinger, H. M. Hofmann, and G. M. Hale, Phys. 
Rev. C 64, (2001) 044003 

17. H. M. Hofmann and G. M. Hale, Phys. Rev. C 68, 
(2003) 021002; Phys. Rev. C 77, (2008) 044002 

18. S. Quaglioni and P. Navratil, Phys. Rev. Lett. 101, 

(2008) 092501 

19. R. Wiringa, private communication 

20. A. Kievsky et al, J. Phys. G: Nucl. Part. Phys. 35, 
(2008)063101 

21. M. Viviani, S. Rosati, and A. Kievsky, Phys. Rev. Lett. 
81, (1998) 1580 

22. M. Viviani et al, Phys. Rev. Lett. 86, (2001) 3739 

23. R.B. Wiringa, V.G.J. Stoks, and R. Schiavilla, Phys. 
Rev. C 51, (1995)38 

24. M. Viviani et al, Few-Body Syst. 39, (2006) 159 

25. L.E. Marcucci et al, Phys. Rev. C 80, (2009) 034003 

26. M.Viviani et al, Few-Body Syst. 45, (2009) 1 19 

27. D.R. Entem and R. Machleidt, Phys. Rev. C 68, (2003) 
041001 

28. E. Epelbaum et al, Phys. Rev. C 66, (2002) 064001 

29. P. Navratil, Few-Body Syst. 41, (2007) 117 

30. V. Bernard et al, Phys. Rev. C 77, (2008) 064004 

31. B.S. Pudliner et al, Phys. Rev. C 56, (1997) 1720 

32. A. C. Fonseca, Phys. Rev. Lett. 83, (1999) 4021 

33. T. W. Phillips, B. L. Berman, and J. D. Seagrave, Phys. 
Rev. C 22, (1980)384 

34. K. F. Famularo et al, Phys. Rev. 93, (1954) 928 

35. D. G. McDonald, W. Haberli, and L. W. Morrow, 
Phys. Rev. 133, (1964) Bl 178 

36. B. M. Fisher et al, Phys. Rev. C 74, (2006) 034001 

37. M. T. Alley and L. D. Knutson, Phys. Rev. C 48, 
(1993) 1890 

38. Y. Koike and J. Haidenbauer, Nucl. Phys. A463, 
(1987) 365c 

39. H. Witala, W. Glockle, and T. Cornelius, Nucl. Phys. 
A491, (1988) 157 

40. A. Kievsky et al, Nucl. Phys. A 607, (1996) 402 

41. T.V Daniels et aL larXiv: 1883 ■ 5 868 

42. T.V. Daniels, private communication 



19 International IUPAP Conference on Few-Body Problems in Physics 



43. F. Zernike and H.C. Brinkman, Proc. Kon. Ned. Acad. 
Wensch. 33, (1935) 3 

44. M. Fabre de la Ripelle, Ann. Phys. (N.Y.) 147, (1983) 
281 

45. M. Abramowitz and I. Stegun, Handbook of Math- 
ematical Functions (Dover Publications, Inc., New 
York, 1970) 

46. M. Viviani, Few-Body Syst. 25, (1998) 177 

47. S. Shimizu et ai, Phys. Rev. C 52, (1995) 1 193 



